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We assess the extent to which the NMSSM can allow for light dark matter in the 2 GeV <, m~o <J 
12 GeV mass range with correct relic density and large spin-independent direct-detection cross 
section, asi, in the range suggested by CoGeNT and DAMA. For standard assumptions regarding 
nucleon s-quark content and cosmological relic density, p, we find that the NMSSM falls short by 
a factor of about 10 to 15 (3 to 5) without (with) significant violation of the current (g — 2) M 
constraints. 
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The CoGeNT collaboration has announced detection of very low energy events which are not consistent with any 
known backgrounds [T] . One possible interpretation of these events is elastic scattering of a light dark matter particle 
(to ~ 5 — 10 GeV) with a spin-independent cross section, a si, on the order of 2 x 10 -40 cm 2 (i.e. 2 x 10~ 4 pb) [T|- 
Hj. This is not very far from the region required to explain the annual modulation observed by the DAMA/LIBRA 
collaboration [3] . A consistent interpretation of both the DAMA and CoGeNT observations [5] is for dark matter to 
Qh have mass and cross section in a 2a ellipse ranging from 0-5/ ~ 3x 10 -4 pb at to ~ 6 GeV down to a si ~ 1.4 x 10 -4 pb 
at to ~ 9 GeV. Clearly, it is of great interest to explore different kinds of dark matter models with regard to their 
qj ability to yield large a si for to ~ 6 — 9 GeV. 

A number of groups have addressed this issue within the context of the minimal supersymmetric standard model 
(MSSM) |H|. However, given the structure of the MSSM Higgs sector and constraints thereon from LEP and 
elsewhere, achieving the above cross section at low LSP mass is not possible [HI IS] • A much higher local density of 
dark matter than the measured cosmological dark matter density, p = 0.3 GeV/cm 3 , would be needed to bring the 
a si required to describe the CoGeNT/DAMA events down to the level possible within the MSSM. Basically, the 
problem is that the Higgs with Enhanced coupling to down quarks, whose exchange is primarily responsible for the 
elastic scattering of the LSP (the lightest neutralino) on a nucleon, must be rather heavy in the MSSM context after 
imposing LEP constraints. Of course, a local density much larger than the cosmological average could be assumed so 
Q\ as to get the needed asi at low to. However, there is a second problem. For low LSP mass, the MSSM simply does 
not allow sufficient early universe annihilation to yield the observed cosmological average relic density once Tevatron 
limits on B(B S —> p + pr) are imposed |S]. 

Thus, it is interesting to see if an extension of the MSSM could allow the relevant Higgs boson to have lower 
mass than allowed in the MSSM, thereby achieving asi — (1.4 — 3.5) x 10~ 4 pb, while maintaining consistency with 
all constraints. In a previous paper [TU], we explored this question within the context of supersymmetric models 
with an additional generic chiral singlet superfield and found that this was indeed possible, the successful scenarios 
being ones in which both the LSP and exchanged Higgs are substantially singlet in nature. In this paper, we focus 
on the concrete (and more restrictive) case of the next-to-minimal supersymmetric standard model (NMSSM) . Our 
conclusion will be that the observed cosmological relic density can be achieved while maintaining consistency with 
limits on B(B S —> p + pr) but that the largest asi values that can be achieved for standard inputs regarding the 
s-quark content of the nucleon fall short of the preferred asi region of [5] by a significant factor. In particular, in 
the strict NMSSM, scenarios with a light singlet Xi an d largely singlet light Higgs cannot be realized at high tan/3 
while satisfying all other constraints. We also briefly discuss possibilities for enhancing the NMSSM cross sections by 
enhancing the s-quark nucleon content or reducing the required asi using the recently proposed larger local density 
p ~ [0.4 — 0.485] GeV/cm 3 (see [TT] for a summary). 

The remainder of this paper is structured as follows. In Sec. |Hj we outline the problems faced in the MSSM. In 
Sec. |III| we discuss how the NMSSM can potentially avoid these problems without violating the relevant collider 
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constraints. In Sec. IV we turn to a detailed discussion of the NMSSM, including the point searching procedures 
we will employ and the constraints that must be obeyed. In Sec. [V] we present the NMSSM benchmark points we 
have found with large asi that satisfy all LEP and BaBar limits. We then examine implications of various additional 
constraints from the Tevatron, B physics and (g — 2) M for such points. We discuss some phenomenological issues for 
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those points that survive all constraints. In Sec. VI we summarize our results and draw conclusions 



II. LIGHT NEUTRALINOS IN THE MSSM 

In the MSSM, there are two CP-even Higgs bosons, the h° and the H° with m^o < m^o . In the usual convention, 
one writes H° = cos aHd + sin aH u , h° = — sin aHd + cos aH u , where Hd, u are the neutral Higgs fields that couple to 
down and up type quarks respectively. An especially crucial parameter of the model is tan/3 = (H u )/(Hj). Relative 
to the SM Higgs, g^vv = sin(/3 — a) and gn°vv = cos(/3 — a), where VV = W + W~~ or ZZ. The structure of the 
model combined with LEP constraints require that m^o, rn^o > 90 — 100 GeV. In this case, cos(/3 — a) must be fairly 
small, especially at large tan/3. The combination of large tan/3 and small cos(/3 — a) implies a ~ and cos a ~ 1. 
In this situation, the only way to get a large spin-independent cross section for lightest neutralino, x?, scattering on 
the nucleon is via exchange of the H between the \i (9h°x°x° ^ COSQ and the down type quarks contained in the 
nucleon (g_H°dd,ss,bb tan/3cosa). A rough formula for the spin-independent cross section was obtained in [TU] : 

,n , n -B v. f N h\ Aan/3\ 2 /100GeV\ 4 4 

* ai w L7xl ° p H^fJl^rJ l^rJ cosV (1) 

where we have written xl = NuB + Ni 2 W 3 + N 13 Hd + Ni 4 H u . In the above, Nf 3 cannot be much larger than 0.1 
because of limits on the Z invisible width. Given that LEP constraints basically force mjjo > 100 GeV and that 
other constraints (including 6-quark Yukawa perturbativity) are very difficult to satisfy for tan/3 > 50, we see that 
the MSSM is unable to obey all constraints and yield &si larger than a fraction of 10~ 4 pb. 

In addition, one must consider whether the MSSM allows for sufficient early- universe annihilation to achieve O^o/i 2 < 
0.1. To briefly review, the density of neutralino dark matter in the universe today can be determined by the particle's 
annihilation cross section and mass. In the mass range we are considering here, the dominant annihilation channel 
is to bb (or to a lesser extent to t + t~) through the s-channel exchange of the pseudoscalar Higgs boson, A. The 
thermally averaged cross sections for these processes are given by 



,2 



_ (3, ljggmg T tan 2 /3 ™gj V 1 ~ m hl m \\ 

x [(N 13 sin p-N 14 cos /3)(g 2 N 12 -g 1 N 11 )} 2 7 (2) 

where T^o is the width of the pseudoscalar MSSM Higgs. And although there are additional contributions from 
scalar Higgs exchange, these are suppressed by the square of the relative velocity of the neutralinos, and thus are 
substantially suppressed in the process of thermal freeze-out. 
The thermal relic abundance of neutralinos is given by 

10 9 my* l 

ft^o h 2 « — ^— -, — (3) 

Xl Mpi \ a 'x 1 x 1 v ' 

where is the number of relativistic degrees of freedom available at freeze-out and Tpo is the temperature at which 
freeze-out occurs: 



Tfo VV 8 ^3 m o/ Tfo 



For the range of masses considered here, and for cross sections which will yield approximately the measured dark 
matter abundance, we find m^o/Tpo ~ 20. 

For m^o ^5—15 GeV, the relic abundance of MSSM neutralinos is then approximately given by 

o u2 mf 0A \f 50 ^V m A° \V9GeV\ 2 

^ h ~ 0A {Wj{^) iMGevj U^J ■ (5) 

Given that LEP limits require m^o > 90 — 100 GeV and that tan j3 as large as 50 is already in the non-perturbative 
domain for the 6-quark coupling, it requires a very extreme choice of parameters to get the measured dark matter 
density of our universe to be as small as that measured, f^cDM^ 2 = 0.1131 ± 0.0042 [15]. And, even with such 
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extreme parameter choices, a si can be no larger than ~ 1.7 x 10 5 pb. Of course, it is true that the same extreme 
choice of parameters that minimizes Q^oh 2 , bringing it close to the observed value, at the same time maximizes asi- 
However, there is a further barrier to achieving the minimal f^o/i 2 , maximal a si scenario. In particular, the above 
discussion does not yet include consideration of the Tevatron limits on B(B S — > In [5] (see their Fig. 3b), 

it was found that the MSSM simply cannot give the correct relic density for m^o in the CoGeNT/DAMA region 
once the B(B S — > limit is imposed in addition to the LEP limits. This situation motivates us to consider 

supersymmetric scenarios beyond the MSSM. In the next section, we will demonstrate that in the NMSSM it is 
possible to alleviate both the elastic scattering cross section and relic abundance problems found in the MSSM. 

III. THE NMSSM 

In the NMSSM, one adds exactly one singlet chiral superfield to the MSSM. As is well known, this allows a completely 
natural explanation for the size of the /i term 13J and can reduce electroweak fine-tuning |14j . and potentially catalyze 
electroweak baryogenesis |15j . The NMSSM superpotential is given by 

AM^ + ^S 3 , (6) 

and the associated part of the soft Lagrangian is given by 

XA x SH u H d + IkA k S 3 + H.c. (7) 

The restriction to the forms given above is implemented by invoking a Z3 symmetry to remove all other possible terms. 
In particular, only the dimensionless A and k superpotential terms are allowed. All dimensionful parameters are 
generated by soft-SUSY-breaking. An effective /i value is automatically obtained as /i c ff = X(S). This very attractive 
extension of the MSSM allows for a considerable expansion of the phenomenological possibilities. In particular, the 
singlet superfield leads to five neutralinos, three CP-even Higgs bosons (/ii,2,3) and two CP-odd Higgs bosons (ai^)- 
In general, the neutralino mass eigenstates are mixtures of the MSSM neutralino fields and the singlino field that is 
part of the singlet superfield; the CP-even (odd) Higgs mass eigenstates are similarly mixtures of the CP-even (odd) 
MSSM fields and the CP-even (odd) components of the complex singlet scalar component of the singlet superfield. 

Within the NMSSM, it is very natural for the lightest pseudoscalar Higgs, a%, to have low mass (see [16 ). In 
particular, U(1)r or C/(1)pq symmetries can appear which lead to values of m ai well below the electroweak scale. If 
one is close to either symmetry limit, the a\ will be at least moderately singlet-like (as opposed to being more purely 
MSSM-Higgs-like) and will likely be beyond the reach of current collider constraints. 

That a light a\ in the NMSSM can allow a very light dark matter particle in the CoGeNT mass region with correct 
relic density was established in |17) . This is because the light a\ s-channel annihilation process is typically fairly close 
to being 'on-pole', 2m^o ~ m ai , as opposed to 2to^o <§; m^o for the rather heavy A of the MSSM. However, in 
the scans performed in |17j we did not encounter points with cross sections as large as those needed to describe the 
tentative CoGeNT/DAMA signal. We now describe a strategy for getting the largest possible cross section. 

To enhance the neutralino's elastic scattering cross section, we need a Higgs mass eigenstate that is primarily Hd 
(so that it will have enhanced couplings to down-type quarks at large tan/3) with mass lower than possible for the 
H° of the MSSM. While this is not as easy to arrange in the NMSSM as are low values of the lightest CP-odd Higgs 
mass, it is still possible. The value of the down-type diagonal term of the NMSSM scalar Higgs (squared) mass matrix 
at tree-level is given by 

(J 2 V 2 

w2 H ^2 = —— — 2~S +Mtan/3(A A +Kfi/X), (8) 
1 + tan^ p 

where v is the Standard Model Higgs vacuum expectation value. At large tan /3, in order for this to fall significantly 
below the value of the up- type Higgs entry (which is generally m 2 H u ss (85GeV) 2 ), there must be some cancellation 
between the A\ and K/i/A terms. This cancellation also suppresses the mixing term between up-type and down-type 
scalar Higgs bosons. The down-type mass can further be protected from large radiative corrections if the two stop 
masses are similar. Together, these features can potentially lead to a down-type scalar Higgs boson with a mass 
significantly below 100 GeV. 

The scenarios that can potentially lead to large a si are then ones in which the lightest of the NMSSM Higgs bosons, 
the hi, is not SM-like, has enhanced coupling to down-type quarks and has mass below ~ 100 GeV. The hi will 
typically be SM-like and for below 100 GeV is typically not very heavy - m^ 2 > 110 GeV for tosusy = 500 GeV 
and rrih 2 115 GeV for msusY = 1 TcV. LEP limits will be very constraining in this situation. In addition, 
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many P-physics constraints will enter as will constraints from (g — 2)^. Also important will be limits on bb + Higgs 
production with Higgs — > t+t~ and t — > h + b decays with h+ t+v t . We will employ augmented versions of 
NMHDECAY [TSJ [TH] supplemented by micrOMEGAs gU] (the latter will be implemented as in NMSSMTools [21]) 
for our exploration of the NMSSM parameter space. 

IV. CONSTRAINTS AND SCANNING IN THE NMSSM 

As noted, we have performed our scanning using an augmented version of NMHDECAY linked to micrOMEGAs 
as in NMSSMTools. NMHDECAY currently incorporates all LEP limits on Higgs bosons as well as LEP limits on 
neutralinos and charginos. 1 We have augmented NMHDECAY to include the recent ALEPH constraints [23] on 
e + e~ — » Z + Higgs with Higgs — > aa (in our case a = ai, the lightest CP-odd Higgs boson of the NMSSM) with 
a — > t+t~ . Further, we have augmented NMHDECAY to include the combined CDF+DO Tevatron constraints [23] 
on bb + Higgs production with Higgs — > t + t~ (for the scans performed in this paper, it is constraints in the case 
of Higgs = hi or a 2 that are typically relevant). 2 Finally, in the scenarios with large a si the h+ is inevitably light 
enough that t — > h + b decays will be present and, since tan/3 is large, h+ — > t+v t will be completely dominant. We have 
thus augmented NMHDECAY to include the current DO limits [35] on B{t ->■ h+b) x B(h+ -> t+v t )? NMHDECAY 
also includes analysis of a large selection of B physics constraints. For our purposes, the most important ones turn 
out to be B s — > Lt+(i~, B+ — » t+v t , and b — > sj. We have also augmented NMHDECAY to incorporate full BaBar 
constraints on T„s — > ja with a — > fi +fj,~~ or a — > t+t~ as implemented in |26j . Finally, we have examined the 
NMHDECAY predictions for (g — 2)^ for high-as/ cases. In our search for desirable points, we have demanded that 
all the LEP limits, including the ALEPH limits, are strictly obeyed. We have also demanded that the BaBar limits 
be strictly satisfied. 

The bb+Higgs(-^ t+t~) and t — > h+{— > T+v)b limits are treated somewhat differently In the experimental papers, 
the observed limits are plotted as a function of the relevant Higgs mass in comparison to the expected limits. The 
expected limits have error bars that are partly statistical and partly systematic (including theory systematics) that 
have been combined in quadrature, i.e. assuming a Gaussian distribution in particular for theoretical systematics. 
We believe that treating the observed limits in these cases as true limits is somewhat dubious. In our opinion, it 
would be much better to have separated the statistical errors from the systematic errors and ask what band about 
the observed limits would result from pushing all systematics in the least or most favorable direction. In the absence 
of sufficient information to carry out this task, we will simply assess the impact of relaxing the observed limits in 
the above channels by an amount equivalent to the \a or 2a error bands (as plotted relative to the expected limits) 
relative to the observed limits. 4 

In assessing the B+ t+v t , b — > sj and (g — 2) M constraints contained in the basic NMHDECAY program 
(B s — > n+n~~ is handled differently as described later) we have adopted the following procedure. The NMHDECAY 
output gives the model point prediction as well as the maximum and minimum values after adding and subtracting 
the theoretical error. Let us call these Po, P + and P_, respectively. Also contained in the output is the ±2<r 
interval for the experimentally observed value or limit, which we label as 0+ 2(T and 0_ 2(T , respectively. Any point 
for which P + or P_ falls within the interval / = [0-2a, 0+2a\ is deemed acceptable. If this is not the case we assess 
the extent of the violation of the constraint as follows. Let us say P_ > 0+2 CT . Define A = |P_ — 0+2a\- We 
then compute R a — A/E, where E is a combined error associated with the experimental and theoretical errors: 
E = [{\0 +2 a - 0_ 2ff |/4) 2 + (|P+ - P-l/2) 2 ] 1 / 2 . If P+ or P_ falls within the interval I = [0_ 2( t, +2ct ] we set R a = 0. 
We will summarize the values found for R a for high-erg/ points for each of the above three constraints. 

In our scans, we have held fixed the soft scales M 2 = 200 GeV and M3 = 300 GeV, allowing for varying values of 
Mi (which essentially fixes the mass of the bino-like neutralino). Our scans have been performed for fixed values of 
/Xes = +200 GeV and —200 GeV. (It seems that smaller |/z ef j| values do not allow large asi to be consistent with 
all other constraints. Conversely, larger |// e ff| tends to lower the achievable a si-) We have considered three values of 
tan f3, tan j3 = 40, tan j3 = 45 (only for /i e g < 0) and tan /? = 50. We have adopted a universal value of msuSY for all 
the soft SUSY-scale slepton and squark SUSY-breaking masses. We consider tosusy = 500 GeV and 1 TeV. We have 
adopted a universal value for all the soft A parameters, i.e. A so m = A t = = A T , .... It turns out that essentially 



1 We have retained the stronger cross section constraints of the original NMHDECAY program rather than weakening them in the manner 
suggested in |22| . However, we have updated the limit on r z _ t -o-o to 1.9 MeV as in 1221 . 

2 Experimental plots assume the MSSM for which the H and A are nearly degenerate whereas in most NMSSM cases /i2 and 0,2 are not 
degenerate, implying a somewhat weaker constraint on the separate bbh,2 and bbai couplings. 

s Limits in this channel from CDF are not currently available. 

4 In the t — ¥ h + b case, plots only show a la error band. We have simply doubled this for an approximation to the 2a error band. 
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the only way to obtain a value for B(B S —> fi + n~) below the current experimental limit when tan/3 is large is to 
choose A so f t rather precisely (typically to within 1%). At high tan/3, it turns out that the appropriate choice for 
A so ft is essentially only a function of msuSY- For each choice of tosusy, we- have determined the appropriate A so f t 
and have then held it fixed at this value as we scan over other parameters and assess all the other constraints (LEP, 
BaBar, Tevatron, . . . ). 

In all our scans, we have consistently found that large asi is only achieved if the Xi is mostly bino, implying that 
m^o is pretty much fixed to be close to Mi. As a result, we have performed scans at a variety of M\ values in the 
general CoGeNT range. For any given M\ we thus end up scanning in A, k, A\, A k , demanding, as sketched above, 
complete consistency with all LEP and BaBar limits, but allowing for some deviation from i?-physics, Tevatron and 
(g — 2) p nominal constraints. For a choice of A, k, A\, A k that is allowed by LEP and BaBar constraints (at the given 
A so ft) i there is no guarantee that Vlti 1 ~ 0.11 will be obtained. Fortunately, it is often the case that one can adjust 
m ai (by changing A K by a relatively small amount) and or m^o (by changing M\) so that flh 2 ~ 0.11 (we accept 
points within the NMSSMTools-defined window, 0.094 < ilh 2 < 0.136) is achieved without destroying consistency 
with LEP and BaBar limits. The results of these scans after this adjustment are presented in the following section. 

V. BENCHMARK MODELS IN THE NMSSM 

We begin with plots, Figs, [l] and [5J of asi vs - m%o for // e g = —200 GeV and [i e g = +200 GeV. We only give points 
found that have fairly large asi- For these two figures, only the LEP constraints, BaBar constraints, B(B S — > /i + /i~) 
limits and 0.094 < fi/i 2 < 0.136 are required to be satisfied. We refer to these as level-I constraints. Many of the 
plotted points with the largest asi values fail at some level one or more of the other constraints, as we shall describe. 

For /i G ff = —200 GeV we see in Fig. [I] that fairly large values of asi (only a factor of 3 to 5 or so below the values 
typical of the preferred CoGeNT/DAMA region) can be obtained. Such points typically have both large tan/3 = 50 
and low msuSY (so that rrih 1 can be relatively smaller). In contrast, Fig. [2] shows that for /j, c g = +200 GeV we 
never found any points with tan/3 = 50 and tosusy = 500 GeV that were consistent with LEP and BaBar limits. 
Consistent points were found for tan/3 = 40 and tosusy = 500 GeV with asi ~ 0.1 x 10~ 4 pb. For msuSY = 1 TeV, 
consistent points are found for both tan /3 = 50 and tan /3 = 40 for which the largest cross sections found are of order 
0.2 x 10 -4 pb and 0.15 x 10~ 4 pb, respectively, both of which are significantly below the cross section needed to 
explain CoGeNT/DAMA events. 

As anticipated from our earlier discussions, one finds that almost all the high-o si points for either sign of /z e ff have 
Cv(hi) <C 1 (where Cy{h) = ghvv / 9h SM vv), implying that either h% or h 3 is the SM-like Higgs boson. This was 
not imposed, but simply came out of the scan when large asi was required. This shows that our intuition as to how 
to achieve large asi was correct. For many cases, m h . 2 < 110 GeV and Cy{h 2 ) ~ 1. Such points escape LEP limits 
because B(h 2 — > a\ai) is large and 10 GeV < m ai < 2ms, the 10 GeV lower bound so that BaBar constraints on 
T3S — > J0>i and ALEPH constraints on Zh 2 with h 2 — >■ a\a\ — > 4t are obeyed and the upper bound so that a\ — > bb 
is forbidden. 

Of interest for the following are the masses of the h 2 and h + for the large a si points. These are shown in Figs. [3] 
and El One should take note of the rather low values of , and m h + . (For some points, to/j 3 is also quite 
smalL) Low combined with large tan/3 implies that B(t — > h + b) will be significant and that B(h + — > t + v t ) ~ 1. 
Low masses for the neutral Higgs bosons coupled with the fact that at least several of them will have enhanced bb 
coupling when tan/3 is large implies that bb + Higgs with Higgs — > t + t~ will have a high rate at a hadron collider 
for several of the neutral Higgs. Thus, Tevatron constraints will often be of importance, and future LHC results could 
have a deciding impact. 

Indeed, let us now add to the LEP and BaBar constraints the requirement that the Tevatron constraints on 
bb + Higgs and t — > h + b be satisfied within la as defined in the previous section. The plot of Fig. [H] shows that for 
/i c ff = —200 GeV the points with largest asi (i-e. those with low msuSY and hence lower and large tan/?) do not 
satisfy the additional Tevatron constraints. The maximal cross section allowed is ~ 0.3 x 10~ 4 pb, which is distinctly 
below the a S i = (1.4 - 3.5) x 10~ 4 pb of the CoGeNT/DAMA region. 

For fi e g = +200 GeV, the LEP and BaBar constraints had already eliminated such points and imposing the Tevatron 
constraints at the la level eliminates only the single point of Fig. pjwith m^o ~ 2.4 GeV and 0-5/ ~ 0.28 x 10~ 4 pb. 

If we require that the Tevatron observed limits apply with no allowance for error, we obtain the plots shown in 
Fig. |6j The maximum 0-5/ for both fi e g = —200 GeV and fi e g — +200 GeV in the CoGeNT m^o region is of order 

0.14 x 10~ 4 pb, a factor of 10 - 20 below the a S i = (1.4 - 3.5) x 10~ 4 pb CoGeNT/DAMA region. 

We now turn to the impact on these results of B + — > t + v t , b — !• s-f and (g — 2) M constraints. To assess these 
impacts, we employ R a defined earlier, where R a is computed for each of the above three cases. For /z e ff = —200 GeV, 
non-zero values of R a only arise for B + t + v t and (g — 2)^. R (T (B + — s- t + v t ) for the plotted points is typically 
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NMSSM Cogent-like points: /x=-200 GeV 
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for fi e g = —200 GeV. Parameters not shown are fixed as stated in the text. Only level-I constraints are 



below, often well below, 0.4, which we do not regard as a significant exception to the experimental limits. On the 
other hand R a {{g — 2) M ) is often quite large. Indeed, if we require R a {{g — 2) M ) < 2 then all points are eliminated 
except for those with very low m^o ~ 2.4 GeV. Requiring R a {{g ~ 2) M ) < 3 leaves the points plotted in Fig. [7J i.e. it 
is the tosusy = 1000 GeV points that can survive this very loose constraint. In short, if (g — 2) M is taken seriously, 
the fi c ff = —200 GeV points must be eliminated from consideration. Of course, one should never completely rule 
out the possibility that significant additional new physics could contribute to (g — 2) M without affecting the NMSSM 
structure of the Higgs and dark matter sectors. 

In contrast, the vast majority of the fi c s = +200 GeV points (and indeed all of those near the CoGeNT mass 
window) are fully consistent with both B + t + v t and (g — 2)^ constraints within the NMHDECAY windows and 
only have small values of R a (b — > S7). For all of the plotted points in the m^o > 4 GeV region, R a (b — > S7) € [0.5, 0.6]. 
Given the possibility of other new physics that might enter into b — > S7 that might easily have no affect on the NMSSM 
Higgs and dark matter issues, we regard this as acceptable. 

Let us focus on a few more details regarding the [i c s = +200 GeV points. As already noted, only these are fully 
consistent with (g — 2) M constraints. As described above, they have only a small violation of nominal b —¥ S7 bounds. 
In the left-hand plot of Fig. |8j we show the range of m ai values as a function of m^o. One observes the expected 
trend of increasing m ai with increasing m^o needed in order to achieve appropriate relic abundance. 

As discussed in Ref. [37], scenarios with a light a± can potentially be probed by directly searching for the a\ at 
hadron colliders. The discovery potential is basically a function of the strength of the a\bb reduced coupling, C aib j J . In 
the NMSSM context, C ai6 j = cos 9a tan j3, where cosOa specifies the amount of the a\ that resides in the MSSM-likc 
doublet sector as opposed to the singlet component: 



ai = cos 6 acimssm + sin 9 A as 



(9) 



In the absence of cos 9a suppression, the a\ would be strongly coupled to down-type quarks proportionally to tan /3. 
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+200 GeV. Parameters not shown are fixed as stated in the text. Only level-I constraints are 



However, many of the points with large a si have cos 9 a values significantly below unity. The right-hand plot of Fig. [8] 
shows |C 6 j| vs. m ai for all the fj, e g — +200 GeV points. We see significant variation of |C 6 j|, but find many points 
with fairly large values. Of course, the larger \C a b g| is, the easier it will be to detect the a± directly in hadronic 



Some of the 



C ai bb\ ValueS 



collisions, for example via bbai production followed by a,\ — > t + t~ or gg 
are sufficiently large that early detection at the LHC might be feasible. 

Another interesting question is how well the points plotted agree with precision electroweak constraints. This can 
be assessed by computing the effective precision electroweak mass defined by 



lnm e //= ^2 \ c v(hi)\ 2 \nm hz 



(10) 



i=l,2,3 



One finds that all fi c g = +200 GeV points with tosusy = 1 TeV have m e ff E [114 GeV, 116 GeV]. In comparison, 
the tosusy = 500 GeV points can have m e ff as low as 100 GeV, thereby achieving excellent agreement with precision 
electroweak measurements. Such points are closely related to the "ideal" Higgs scenarios, but are more complex in 
nature. However, as shown in Fig. [2j the largest a si that can be achieved for such points is of order 0.1 x 10~ 4 pb, a 
factor of > 15 below that needed to most naturally describe the CoGeNT/DAMA observations. 

Although we have not explicitly performed the necessary computations, we anticipate that the to-susy = 1000 GeV 
points will have significant electroweak symmetry breaking (EWSB) finetuning (i.e. to predict the correct value of mz 
will require very precise adjustment of the GUT-scale soft-SUSY-breaking parameters) whereas much less finetuning 
should be required in the case of the Tosusy = 500 GeV points. 

It is perhaps interesting to give details for the tan/3 = 40, tosusy = 500 GeV "semi-ideal-Higgs" point with m^o in 
the center of the CoGeNT mass region and asi ~ 0.1 x 10~ 4 pb found in Fig. [2} The relevant details are presented in 
Table |lj For this point it is the /12 with m^ 2 ~ 97 GeV that is mainly responsible for a substantial size for a si (since 
^h 2 bb ^ S l ar § e — see the 3rd row of Table H. In contrast, the h% has relatively small down-type quark coupling as can 
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constraints are imposed. There is a lot of point overlap in this plot. 
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be seen from the tabulated value of C yi ^. Note that low m e ff is achieved despite the fact that the Higgs, namely 
the /13, that carries the bulk (74%) of the WW, ZZ coupling-squared has mass m,h 3 ~ 126 GeV. This is because the 
hi carries about 25% of the WW, ZZ coupling-squared and has very low mass. 
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According to the NMSSMTools package, the only statistically significant Higgs signal for this point in the normal 
LHC search channels arises in the WW —> hs —> t + t~ channel where one finds statistical significances relative to 
background of 3.8 and 14 at low and high luminosity, respectively. Even though the h 3 in this scenario is fairly SM-likc 
{\Cv{hz)\ 2 ~ 0.72) its decays to WW,ZZ and 77 are suppressed to levels well below those typical of the SM Higgs 
of the same (low) mass, partly because of the smaller |Cy(/i3)| 2 and partly because of significant h% — > Higgs pair 
decays. In addition to the WW — > h 3 — > t + t~ LHC signal, it seems to us that the bbh 2 (^ t+t~) signal would also 
be strong. One would also wish to push discovery of the a\ in the gg — > ai — > \i + [i~ channel — the preliminary 
estimates of [37] indicate this signal might well be observable given the relatively large value of C ^ tabulated above, 
despite the fact that m ai is in the Upsilon mass region. 

This and other similar points for which /13 is the SM-like Higgs appear distinctly in Fig. [9] In these scenarios, LEP 
constraints are easily evaded for the h\ and h 2 since they have greatly reduced WW, ZZ coupling, and in the case of 
the hi the dominance of hi — > a%ai — > 4r decays greatly reduces LEP sensitivity as well. LEP constraints on the 
do not enter since to^ 3 > 114 GeV for these cases. 

As regards the B physics results in the last row of Table |TJ the possible range of predictions is that obtained by 
taking the central prediction of the point after subtracting or adding the theoretical error. These ranges can be 
compared to the current ±2cr experimental ranges of Table [il] For all but B(b —> 57) there is satisfactory overlap of 
the predicted range with the experimental range. If we quantify the discrepancy between the predicted and observed 
ranges as described earlier, the overlap failure is at about the 0.5cr level. 

Let us briefly discuss the spin-dependent cross sections for the /i = +200 GeV points. These are basically only 
a function of tan/3 and m^o. The proton and neutron spin-dependent cross sections are very similar in magnitude. 
Thus, we confine ourselves to plotting the average value asD = [<^sd + <t s_d)/2 in Fig. 10 even though it is only 
the separate cross sections that are directly experimentally measurable. One finds that a$D varies from a low near 
0.24 x 10 -4 pb for m^o ~ 2.5 GeV to a high of ~ 0.6 x 10~ 4 for m^a ~ 11 GeV (the largest value we have considered 
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for n c ff = +200 GeV points). 

All the cross section results obtained above are based on the nominal NMSSMTools and micrOMEGAs assumptions. 
It is worth mentioning several means of enhancing these cross sections. First, we note that the cross section magnitudes 
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TABLE I: Properties of a particularly attractive but phenomenologically complex NMSSM point with /i = +200 GeV, tan /3 = 40 
and msuSY = 500 GeV. This point predicts values for B(t — > h + b) x B(h + — > t + l> t ) and for bb + Higgs production with 
Higgs — > t + t~ (for all neutral Higgs bosons) below current observed Tevatron limits. In the last row, the brackets give the 
range of predictions for this point after including theoretical errors as employed in NMHDECAY. 
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TABLE II: The ±2<r experimental ranges for the B physics observables tabulated in the last row of Table [I] 
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have assumed the standard s-quark content for the proton. In |22j . the possibility of enhancing 175/ by increasing 
the s-quark content of the nucleon was discussed. In particular, if one changes the nominal micrOMEGAs values of 
(t^^y = 55 MeV, (T = 35 MeV to cr^jv = 73 MeV, cr = 30 MeV then o 5/ will be enhanced by roughly a factor of 
3.3. We believe that such a large shift is not consistent with current constraints and lattice calculations. At most, 
one might consider ovyy ~ 60 MeV and ctq = 30 MeV [28], leading to an enhancement of about 50%. In fact, the 
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FIG. 9: |CV(3)| vs. |CV(2)| and m,h 3 for p c ff — +200 points. Parameters not shown are fixed as stated in the text. Only level-I 
constraints are imposed. There is a great amount of point overlap in this plot. 



preponderance of information suggests that, if anything, a lower value of o^n ~ 50 MeV is preferred leading to a 
decrease in the nucleoli's s-quark content and thereby a decrease in a si- Another possibility is to employ the larger 
average local dark matter density p = [0.4 — 0.485] GeV/cm 3 suggested in recent papers (see the summary of [TT] ) 
instead of the micrOMEGAs default value of 0.3 GeV/cm 3 . This would result in a ~ 60% decrease in the asi required 
to explain the CoGeNT/DAMA events. Using both a 50% s-quark enhancement and the larger p one could get about 
a factor of 2 decrease in the discrepancy between the NMSSM predictions for asi an d the asi values needed to 
describe the CoGeNT/DAMA observations. 

For nominal s-quark content, our results differ somewhat from the NMSSM scan performed in [22] . Their results for 
asi for p e g > are roughly a factor of 10 below ours. We believe that this is primarily because in their scenarios the hi 
is always SM-like, whereas in our highest-os/ cases the hi has enhanced down-type quark coupling and it is the hi or hj, 
that is SM-like. This means that their non-SM-like mainly if^-likc Higgs, the hi in their case, is typically significantly 
heavier than in our scenarios. Since a si oc l/ m fr -Mfte a f ac tor of 10 increase in asi can be achieved if mn d -uke is 
decreased by a factor of 1.77. For their scans, the largest asi is achieved for m Hd _i ike G [205 GeV, 260 GeV], whereas 
our large asi values typically have mH d -uke < 100 GeV. They obtain some gain in cross section since their typical 
Peg is lower (~ 138 GeV vs. our 200 GeV), leading to somewhat larger Nf 3 . The larger mti d -iike in the scans of [22] 
imply a larger mh+ with the consequence that their largest asi points are within the nominal ±2<r constraints from 
b — > S7 whereas our high-<7s/, p c g > points are about 0.5cr outside the ±2cr region. 



VI. SUMMARY AND CONCLUSIONS 



We have examined parameter choices within the NMSSM that are potentially capable of yielding a large spin- 
independent cross section for nucleon-LSP scattering at low LSP mass, consistent with that needed to describe the 
CoGeNT/DAMA observations. We have required that all LEP and BaBar constraints be satisfied and that accepted 
points have correct relic density and sufficiently small B(B S — > p + p~). We have then examined the impact of 
additional constraints associated with Tevatron observations, other B physics observations and (<? — 2)^. 

For standard assumptions regarding the s-quark content of the nucleons, we have found that in the NMSSM the 
largest spin-independent cross section that can be achieved for a relevant range of to^o if p c g > is roughly a factor 
of 10 to 20 shy of that needed to describe the CoGeNT/DAMA event excesses assuming standard relic density, the 
latter corresponding to asi ~ (1.4 — 3.5) x 10~ 4 pb. In particular, asi for p c g = +200 GeV can be no larger 
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than 0.14 x 10~ 4 pb after imposing the Tevatron constraints (but allowing for a very mild violation in b — > S7). If 
one allows for p ~ [0.4 — 0.485] GeV/cm 3 instead of 0.3 GeV/cm 3 this will decrease the asi required to explain 
CoGeNT/DAMA by about 60% to perhaps as low as ~ 10~ 4 pb. Nonetheless, our maximal asi values, of order 
0.14 x 10~ 4 for n e s = +200 GeV, would still be well shy of that needed. There is also some uncertainty in the s-quark 
nucleon content. It is possible to suppose that it could be enhanced by about 50%, although a 50% decrease is perhaps 
even more reasonable. Combining a 50% increase with the larger p, one would still be a factor of at least 5 short of 
explaining the CoGeNT/DAMA event rates. 

For standard s-quark nucleon content, the largest 175/ values found for /x c ff < are ~ 0.6 x 10~ 4 pb, within a factor 
of 3 to 5 of the needed (assuming nominal p = 0.3 GeV/cm 3 ) asi = (1.4 — 3.5) x 10~ 4 pb. Unfortunately, p, c g < 
NMSSM parameter choices yielding such large 175/ all predict an anomalous magnetic moment for the muon that is 
strongly discrepant with the observed (g — 2) M . Nonetheless, it is not impossible that there is some resolution of this 
disagreement coming from physics beyond the NMSSM. 

We have illustrated that Tevatron (and, presumably soon, the LHC) constraints on bb + Higgs production and 
t — > h + b decays are highly relevant for constraining large-erg/ scenarios. Thus, it is clear that if the CoGeNT 
observations really are dark matter detection and if the NMSSM is the relevant model, detection of one or more of the 
ai, hi, ai and h + of the NMSSM at the Tevatron and LHC should be close at hand in the above channels. However, 
it is also the case that detecting the SM-likc Higgs of these scenarios will be very difficult. 

On another front, in a companion paper ]T0] we have demonstrated that allowing an extension of the NMSSM 
to include additional superpotential terms and/or soft-SUSY-breaking terms (while still keeping just one singlet 
superfield) will be sufficiently less constraining that asi values large enough to describe the CoGeNT excess can be 
achieved without any LEP, Tevatron, BaBar, B-physics (other than a quite small b — > sj deviation) or (g — 2)^ issues, 
and using nominal s-quark nucleon content and standard relic density p. The key new feature is that the additional 
parameters allow scenarios consistent with all constraints for which the \i i s highly singlet and the hi is largely 
singlet-like with large XiXi^i coupling and low m/^. 
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